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The equation of state inside very compact objects like neutron stars is still largely unkown. Even 
though a lot progress has been made in recent years to develop the so-called realistic equations of 
state, a lot of insight can be gained by using polytropic equations of state to integrate the stellar 
equations of structure. In this work we provide a brief review of the Newtonian and relativistic 
equations of structure and present some numerical results, which we believe that can be useful 
for students starting to work on this field. The internal structure of the Newtonian polytropes 
is obtained by the numerical integration of the Lane-Emden equation, and we used our results 
which can also be used to study the Chandrasekhar limit for white dwarfs. However, Newtonian 
physics cannot correctly describe very compact and massive objects. There is an upper mass limit 
for neutron stars, supported by observations, which is not predicted by the Newtonian equations. 
Neutron stars are best described under the framework of General Relativity. The introduction of 
the TOV equations (as well as the relativistic Lane-Emden equations) is, therefore, necessary to 
correctly identify stable and non-stable models via the mass-radius relation. The analysis of the 
EoS also becomes relevant to discard models that allow a possible violation of causality (sound 
speed larger than c). 
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Stabihty. 
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1. Introduction 



The current period of discoveries in the fields of Physics allows the acquisition of knowl- 
edge and understanding of physical phenomena previously unknown. Cosmology, for instance, is 
no different. The accelerated expansion of the universe, the recurrent related terms energy and 
dark matter, the possible emission of gravitational waves in binary systems among others ai^e the 
new phenomenology present nowadays. Amongst this phenomenological diversity, one can also 
highlight an object of fundamental importance to the physical verification of gravitational waves 
emission: neutron stars. Thus, the study of these compact objects is important for the substantiation 
of predictions made based on the theory of General Relativity [|l]]. 

Accurate models of neutron stars are still under development. The equation of state of matter 
in the extreme conditions found in the center of those objects remains unknown. There is a variety 
of factors that can be taken into account. The existence of a solid crust, magnetic fields, rotation 
and superfluidity are just a few of them[|l|]. The study of neutron stars is therefore highly interdis- 
ciplinary, covering many areas of Physics. Nevertheless, it is necessary to check whether models 
derived from different assumptions and simplifications are physically consistent and correct. Due 
to its high compactness it is natural that neutron stars are studied under the framework of General 
Relativity. 

In this work we have studied the non-relativistic (Newtonian) and relativistic stellar equations 
of structure with a polytropic equation of state, and numerical results were obtained with a MAT- 
LAB code. Although these results are not new in the literature, we believe that the compilation of 
the relevant equations and some numerical results will be useful for students starting on this field. 

The paper is organized as follows. Section ^ briefly review the relevant equations for stellar 
structure used in this paper In Section ^ we present the numerical results obtained for relativistic 
and non-relativistic polytropes and finally in Section ^ we present our final remarks. 



2. Stellar equations of structure 

2.1 Non-relativistic stars 

A non-relativistic star in equilibrium can be described by the equations 

dP _ GM{r)p{r) 
dr 

M{r) = f Anpr^dr, (2.2) 
Jo 

where r is the radial coordinate inside the star, P is the density, p is the total mass density and M{r) 
is the total mass inside a radius r. An equation of state P{p) also needs to be supplied, and we use 
a polytropic equation of state of the form 

P = Kp^, 7=1 + - (2.3) 
n 

where n the polytropic index and K may be established by specifying the density and pressure at 
the center 
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Introducing the new variables x and y defined by the system 



X 

r = — 



A' 

Pcf{x), 



(2.4) 
(2.5) 



where 



{n+\)K 



1/2 



one can use eqs. ([2. 1[), ( |2.2| ) and ( |2.3| ) to obtain the Lane-Emden equation: 



(2.6) 



(2.7) 



The Lane-Emden functions are the solutions of eq. ( p7| ) with initial conditions y{0) = 1 and 
g (0) = 0. On the surface of the star x{R) = AR and y{R) = 0. 

After solving the Lane-Emden equation ( p7[ ) for a given value of the polytropic index n, we 
can fix the radius R and the total mass M{R) of the star. The parameter A is fixed by the relation 
x{R) = AR, the central density is given by 



M{R)A^ 



4Kx^{R)y'{R) 



(2.8) 



and the constant K can be obtained from eq. ( g^ . Finally, p(r) and p{r) are given by eqs.(p3|) 
and(0). 



2.2 Relativistic stars 



In the relativistic treatment, eq. (gjj) is replaced by the Tolman-Oppenheimer-VoUcoff (TOV) 
equation given by 



dP GM + 4kGPP/c^ 



P 

p + — 



(2.9) 



dr r'^{l-lGM/c^r) 
The relativistic Lane-Emden equation can be obtained by introducing the new dimensionless 
coordinates (similarly to what was done in the non-relativistic case): 



47rp, 
P_ 

Pc ' 



m{r) , 



(2.10) 
(2.11) 

(2.12) 



With the definitions ( 2.1(]| ), (2.11) and (2.12) and the equation of state (2.3), eqs. (2.2) and 
(p!^) can be rewritten as 



dv _ .2 n 

.2<i0 _ (v + a0<^t/v/J(§)(l + a0) 
' ^ ~ l-2a(« + l)v/,§ 



(2.13) 
(2.14) 
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n 


x{R) 


y'{R) 





2.4502 


-0.8167 


1 


3.1400 


-0.1061 


2 


4.3526 


-0.1273 


3 


6.8966 


-0.0424 


4 


14.9716 


-0.008 



Table 1: Values obtained for x{R) and y'{R) from the integration of the Lane-Emden equation. 



where 



Pc Kp 



xjn 



(2.15) 



PcC 

and Pc is the central pressure. In the non-relativistic limit cr ^ 0, eqs. (2.13) and (2.14) can be 
recombined to give the Lane-Emden equation (2.7) with x = E, and y = d. 

The relativistic Lane-Emden functions are the solutions of eqs. ( 2. 13| ) and (2.14) with the 
initial conditions 0(0) = 1 and v(0) = 0. The surface of the star is represented by the first zero £,\ 
of 0. 



Having solved eqs. ( |2.13| ) and ( [2. 14[ ) for given values of n and a, one can choose a value of 
to fix the other parameters of the solution and the equation of state. 



3. Numerical results 



We developed a MATLAB code to solve the Lane-Emden equation ( |2.7[ ) and obtain the stellar 
structure inside the polytrope |^. Some relevant quantities obtained for the Lane-Emden functions 
with different values of n are listed on table |l]. 

On figure we present some typical results obtained for a « = 3 polytrope (usually refered to 
as the standard model), representing a star of IM© and \Rq, with Pc = 16,46g/cm^ and K = 3.84 x 
10^"^ dyncm^ /g'^^^ . We point out that a white dwarf can also be modelled as a « = 3 polytrope. In 
a white dwarf, the electron degeneracy pressure provides the means for stabilizing the star against 
its gravitational radiation, and the equation of state of a relativistic degenerate electron gas is of the 
from p = Kp^/^. Note that from eqs. (|j8|) and (2^) it can be shown that for « = 3 the total mass 



M{R) is independent from the central density p^, which corresponds to the Chandrasekhar limit 
for white dwarfs These stars have a mass less than the Chandrasekhar limit of approximately 
I.4M0 and radius of the order of lO'* km [Q]. 

After a low or medium mass star has finished burning most of its hydrogen and helium, it 
can become a white dwarf. At this stage, if the star has a mass below that of the Chandrasekhar 
limit, the collapse is limited by the degeneracy pressure of electrons (the Pauli exclusion principle), 
resulting in a stable white dwarf. However, if the star does not have any more fuel to maintain its 
production of energy and has a mass greater than the Chandrasekhar limit, the pressure exerted by 
the electrons is not sufficient to withstand the force of gravity and the star collapses. Its density 
will increase dramatically, leading to the formation of a neutron star or black hole ' 0]. The 



'it is also possible the formation of a quark star. However, this is a theoretical solution and will not be addressed in 
this project. 
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Chandrasekhar limit is a result from the effects predicted by Quantum Mechanics considering the 
behavior of electrons that cause the degeneracy pressure of a white dwarf. Electrons, as fermions, 
cannot occupy the same energy level, i.e., cannot be described by the same quantum numbers. 
Thus, in white dwarfs, a large amount of electrons are at higher energies, causing a certain amount 
of pressure capable of sustaining the star. For a mass greater than the limit, the degenerescence 
pressure becomes insufficient to prevent the imminent contraction of the star. So far, no white 
dwarf with mass greater than the Chandrasekhar limit has been observed 





3 4 5 6 7 

R(cm) .f,^" 



Figure 1: Numerical results obtained for y{x) (solid line) and dy/dx (dotted line) with n = 3, and radial 
behavior of the density p, pressure p and mass M of a star with IMg and IRq. 



In the relativistic treatment, we have extended our calculations, including all the results pre- 
viously done, to stiffer (incompressible matter, n ^ 0) EoS. The results obtained are presented 
in tables For any solution with given values of n, a and and a chosen value of p^, we can 



determine K from eq. (2.15), and we can obtain R and M{R) with 



R 

M 

M 



1/2 



M. 



(3.1) 
(3.2) 

(3.3) 



47iG 

From these results it is possible to verify which solutions are stable and causal. The velocity 

dP 
dp 



of low frequency sound waves is given by = = ^ = c^y{P / pc^). At the stellar center where 
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a 


^1 


v(^i) 


M 


a 


^1 


v(^i) 


M 


0.0 


6.8966 


2.0182 


2.0182 


0.0 


4.3531 


2.411 


0.0 


0.1 


6.8264 


1.0785 


1.0785 


0.1 


3.6989 


1.2987 


0.4107 


0.2 


7.9494 


0.713 


0.713 


0.2 


3.3983 


0.8403 


0.3758 


0.3 


10.8346 


0.5386 


0.5386 


0.3 


3.271 


0.6055 


0.3316 


0.4 


17.8244 


0.4516 


0.4516 


0.4 


3.2473 


0.468 


0.296 


0.5 


37.2163 


0.4214 


0.4214 


0.5 


3.2967 


0.38 


0.2687 


0.6 


91.0674 


0.4493 


0.4493 


0.6 


3.3986 


0.3201 


0.2479 


0.7 


162.6175 


0.5266 


0.5266 


0.7 


3.5469 


0.2773 


0.232 


0.8 


187.233 


0.5969 


0.5969 


0.8 


3.7334 


0.2457 


0.2197 


0.9 


187.029 


0.6375 


0.6375 


0.9 


3.9539 


0.2216 


0.2102 


1.0 


183.6571 


0.657 


0.657 


1.0 


4.2027 


0.203 


0.203 



Table 2: Lane-Emden relativistic functions parameters for « = 3.0 (left) and n — 2 (right). 

P/pc^ achieves its maximum value a, = c^/a. If a > I/7, the speed of sound would exceed 
that of light and the fluid inside the star could violate causality. 

Both M and M increase with a until a value Ocr, as can be seen in figure ^ for n = \. Since 
dm/dpc > is a necessary condition for stability, Gcr marks the onset of the first mode of radial 
instability [g]. 

0.2502 
0.25 
0.2498 

03 

= 0.2496 
S 

0.2494 
0.2492 
0.249 

0.4 0.42 0.44 0.46 0.48 0.5 0.52 

Sigma 

Figure 2: M as a function of a for « = 1. We can see that Gcr ~ 0.42 separates stable (a < Gcr) and 
unstable models (cJ > <Jcr) 




4. Final remarks 

The present work enabled the acquisition of fundamental knowledge to develop models of 
Newtonian and Relativistic stars. Through the implementation of numerical methods along with 
literature review, it was possible to obtain results that simulate the intra-stellar medium and to 
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o 






M 


G 


^1 


v(^i) 


M 


0.0 


3.1416 


3.1416 


0.0 


0.0 


2.7523 


3.786 


0.0 


0.1 


2.5989 


1.7514 


0.1751 


0.1 


2.2898 


2.1735 


0.1222 


0.2 


2.277 


1.1426 


0.2285 


0.2 


2.0008 


1.4358 


0.192 


0.3 


2.0642 


0.8192 


0.2457 


0.3 


1.8013 


1.035 


0.2298 


0.4 


1.9133 


0.6249 


0.25 


0.4 


1.6544 


0.791 


0.2516 


0.5 


1.8008 


0.4981 


0.2491 


0.5 


1.5409 


0.6295 


0.2647 


0.6 


1.7143 


0.4101 


0.246 


0.6 


1.4503 


0.5164 


0.2727 


0.7 


1.6458 


0.3461 


0.2423 


0.7 


1.3759 


0.434 


0.2779 


0.8 


1.5903 


0.2979 


0.2383 


0.8 


1.3136 


0.3717 


0.2812 


0.9 


1.5447 


0.2605 


0.2344 


0.9 


1.2604 


0.3231 


0.2833 


1.0 


1.5066 


0.2307 


0.2307 










3: Lane-Emden relativistic functions parameters for n = 


1.0 (left) and n = 0.5 


a 


^1 


v(^i) 


M 


a 


^1 


v(^i) 


M 


0.0 


2.504 


4.6089 


0.0 


0.0 


2.4508 


4.9071 


0.0 


0.1 


2.1037 


2.7308 


0.09689 


0.1 


2.0642 


2.9318 


0.0927 


0.2 


1.8438 


1.8377 


0.1781 


0.2 


1.8113 


1.9807 


0.1772 


0.3 


1.6593 


1.3367 


0.2333 


0.3 


1.6306 


1.4452 


0.2375 


0.4 


1.5202 


1.0254 


0.2716 


0.4 


1.4938 


1.111 


0.2811 


0.5 


1.4108 


0.8172 


0.2991 


0.5 


1.3856 


0.8867 


0.3135 


0.6 


1.322 


0.671 


0.3199 


0.6 


1.2976 


0.7283 


0.3385 


0.7 


1.2482 


0.5634 


0.3359 


0.7 


1.2241 


0.6114 


0.3581 


0.8 


1.1856 


0.4816 


0.3485 


0.8 


1.1616 


0.5225 


0.3739 


0.9 


1.1317 


0.4176 


0.3585 


0.9 


1.1078 


0.4532 


0.3869 










1.0 


1.0606 


0.3977 


0.3977 



Tkble 4: Lane-Emden relativistic functions parameters for n = 0.1 (left) and n = 0.0 (right). 



verify the main physical properties of polytropic stars. As for the simulations, the results obtained 
show a very good agreement when compared to previous results found in literature. The classic 
study of the polytrope of index n = 3 allowed us to determine important physical characteristics 
such as the maximum mass limit and relate it to a phenomenon well known to the white-dwarfs, 
the Chandrasekhar Umit. From this, it was possible to understand the concept of degenerate matter 
and its relation to gravitational collapse. 

However, the verification that Newtonian equations do not correctly explain the behavior of 
stellar masses for different radii has become essential. Newtonian equations describe the behavior 
of masses through an indefinite growth whereas General Relativity provides a maximum value of 
mass for each equation of state. It is, therefore, necessary to use relativistic corrections for an 
accurate simulation of compact and massive objects. 

The execution of the simulations allowed, therefore, the development of technical and intu- 
itive notions about the general behavior of stars, understanding the importance of the introduction 
of relativistic corrections and the results arising therefrom, which are essential for determining 
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subluminal and stable solutions. 



References 

[1] B. Schutz, Gravity from the ground up, Cambridge, 2003. 

[2] W.H. Press, B.P. Flannery and S.A. Teukolsky, Numerical Recipes in Fortran 77: The Art of Scientific 
Computing, Cambridge, 1986. 

[3] W.J. Maciel, Introdugao a estrutura e evolugao estelar, Edusp, Sao Paulo, 1999. 

[4] R.R. Silbar and S. Reddy, Neutron stars for undergraduates students. Am. J. Phys. 72, 7 (2004). 

[5] B. Ryden, Introduction to Cosmology, Pearson Addison Wesley, San Francisco, 2003 

[6] J.V. NarUkar, An Introduction to Cosmology, Cambridge, 2002. 

[7] R. D'Invemo, Introducing Einstein's relativity. New York: Oxford, 1992. 

[8] K. C. Chung, Introdugao a Fisica Nuclear, Rio de Janeiro: Eduerj, 2001. 

[9] S.A. Bludman, Stability of General-Relativistic Polytropes, Am. J. Phys. 183, 7 (1973). 



8 



